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We report on a search for large extra dimensions in a data sample of approximately 1 fb _1 of 
pp collisions at y/s — 1.96 TeV. We investigate Kaluza-Klein graviton production with a photon 
and missing transverse energy in the final state. At the 95% C.L. we set limits on the fundamental 
mass scale Mo from 884 GeV to 778 GeV for two to eight extra dimensions. 
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Arkani-Hamed, Dimopoulos, and Dvali (ADD) [l( 
made the first attempt to solve the hierarchy problem 
of the standard model (SM) by postulating the exis- 
tence of n new large extra spatial dimensions (LED). 
In this approach, the SM particles are confined to a 3- 
dimensional brane while gravity is diluted in the larger 
volume. The size of the compactified extra space (R), 
the effective Planck scale in the 4-dimensional space-time 
(Mpi), and the fundamental Planck scale in the (4 + Tri- 
dimensional space-time (Md), are related by the expres- 
sion Af|, ; = SnM^ R n . Due to the compactification of 
the extra space, the gravitational field appears as a se- 
ries of quantized energy states, which are referred to as 
Kaluza-Klein modes. A Kaluza-Klein graviton (Gkk) 
behaves like a massive, non-interacting, stable particle 
whose direct production gives an imbalance in the final 
state momentum as its collider signature. 

In this Letter we report the results of a search for LED 
in the final state with a single photon plus missing trans- 
verse energy (7 + $ T ), using data collected with the DO 
detector at the Fermilab Tevatron collider. This signa- 
ture arises from the process qq — ► jGkk, which is studied 
in detail in 0] • The CDF collaboration carried out a sim- 
ilar search with 87 pb~ x of data, setting 95% C.L. lower 
limits on M D of 549, 581, and 601 GeV for 4, 6, and 8 
extra dimensions, respectively Q. Searches for LED in 
other final states have been performed by collaborations 
at the Tevatron HQ and the CERN LEP collider 0. 

The background to the 7 + $ T signal is dominated 
by electroweak boson production and non-collision back- 
ground where muons from the beam halo or cosmic rays 
undergo bremsstrahlung and produce an energetic pho- 
ton. The electroweak background is dominated by the 
processes Z + 7 — > vv + 7, W — > ev where the electron 
is misidentified as a photon, W + 7 where the lepton 
from the W boson decay is not detected, and W/Z + jet 
production where the jet is misidentified as a photon. 

The DO detector Q comprises a central-tracking sys- 
tem with a silicon microstrip tracker (SMT) and a central 
fiber tracker (CFT), both housed within a 2 T supercon- 
ducting solenoidal magnet, with designs optimized for 
tracking and vertexing at \rj\ < 3 and In | < 2.5, respec- 
tively, where rj is the pseudorapidity [8j measured with 
respect to the geometrical center of the detector. The 
central prcshower system (CPS) is located in front of a 
liquid- argon/uranium calorimeter and consists of three 
layers of scintillating strips, providing precise measure- 
ment of electromagnetic (EM) shower positions. The 
calorimeter has a central section (CC) covering \tj\ < 1.1, 
and two end sections (EC) that extend coverage to \r)\ w 
4.2 Each part contains an EM section closest to the 
interaction region followed by fine and coarse hadronic 
sections. The EM section has four longitudinal layers 
and transverse segmentation of 0.1 X 0.1 in 77 — <f> space 
(where <f> is the azimuthal angle), with the exception of 
the third layer, where it is 0.05 x 0.05. Additionally, 



scintillators between the CC and EC cryostats provide 
sampling of developing showers for 1.1 < \<q\ < 1.4. The 
outer muon system, covering \r)\ < 2, consists of a layer 
of tracking detectors and scintillation trigger counters in 
front of 1.8 T iron toroids, followed by two similar layers 
after the toroids. The data in this analysis were recorded 
using triggers requiring at least one energy cluster in the 
EM section of the calorimeter with transverse momentum 
Pt > 20 GeV. The triggers are almost 100% efficient to 
select signal events. This set of data corresponds to an 
integrated luminosity of 1.05 ± 0.06 fb" 1 [ic|. 

We identify a reconstructed calorimeter cluster as a 
photon when it satisfies the following requirements: (i) 
at least 90% of the energy is deposited in the EM section 
of the calorimeter; (ii) the calorimeter isolation variable 
1 = [£ to t(0.4) - £ cm (0.2)]/£ cm (0.2) is less than 0.07, 
where £'tot(0.4) denotes the total energy deposited in the 
calorimeter in a cone of radius 1Z = \J (A7/) 2 + (A(/>) 2 = 
0.4, and £ , e m(0.2) is the EM energy in a cone of radius 
1Z = 0.2; (iii) the track isolation variable, defined as the 
scalar sum of the transverse momenta of all tracks that 
originate from the interaction vertex in an annulus of 
0.05 < 1Z < 0.4 around the cluster, is less than 2 GeV; 
(iv) it has \rj\ < 1.1; (v) both transverse and longitudinal 
shower shapes are consistent with those of a photon; (vi) 
it has neither an associated track in the central track- 
ing system nor a significant density of hits in the SMT 
and CFT systems consistent with the presence of a track 
with pt in agreement with its transverse energy; and (vii) 
there is an energy deposit in the CPS matched to it. Jets 
are reconstructed using the iterative midpoint cone algo- 
rithm [ll1 | with a cone size of 0.5. The missing transverse 
energy is computed from calorimeter cells with \r)\ < 4 
and corrected for the EM and jet energy scales. 

The photon sample is obtained by selecting events with 
only one photon with pt > 90 GeV, at least one recon- 
structed interaction vertex consistent with the measured 
direction of the photon (see below), and $ T > 70 GeV. 
Additionally, in order to avoid large ]ft T due to mismea- 
surement of jet energy, we require no jets with px > 
15 GeV. The reduction of the signal efficiency due to 
the jet veto on initial state radiation has been estimated 
using PYTHIA [l]| to be about 9%. The applied $ T re- 
quirement guarantees negligible multijet background in 
the final candidate sample while being almost fully effi- 
cient for signal selection. 

We reject events with reconstructed muons and with 
cosmic ray muons identified using the timing of the signal 
in the muon scintillation counters or by the presence of 
a characteristic pattern of hits in the muon drift cham- 
bers that is aligned with the reconstructed photon. In 
order to further reject events with leptons that leave a 
distinguishable signature in the tracker but that are not 
reconstructed in the other subsystems of the detector, 
we impose a requirement on the pr of any isolated track 
not to be greater than 6.5 GeV. A track is considered 
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to be isolated if the ratio between the scalar sum of the 
transverse momenta of all tracks that originate from the 
interaction vertex in an annulus of 0.1 < 1Z < 0.4 around 
the track and the pt of the track is less than 0.3. 

The EM pointing algorithm allows calculation of the 
direction of the EM shower based on the transverse 
and longitudinal segmentation of the calorimeter and 
preshower systems. EM pointing is performed indepen- 
dently in the azimuthal and polar planes. The former 
results in the measurement of the distance of closest ap- 
proach (DC A) to the z axis (along the beam line), and the 
latter in the prediction of the z position of the interaction 
vertex in the event, both with a resolution of about 2 cm. 
We require that the z coordinate of at least one interac- 
tion vertex in the event be within 10 cm of the position 
predicted by the pointing algorithm and use the DCA 
to estimate the remaining background from jet-photon 
misidentification and non-collision events. Misidenti- 
fied jets have poor pointing resolution, and therefore a 
wider DCA distribution compared to electrons or pho- 
tons. Likewise, one can anticipate the DCA distribution 
for photon candidates in non-collision events to have an 
even wider shape. After these requirements, 35 events 
are selected in the photon sample. 

We prepare three DCA distribution templates: the 
non-collision template, the misidentified jets template, 
and the e/7 template. The first template is obtained 
from a sample in which a photon candidate, passing the 
same quality requirements as for the photon sample, is 
selected from events with no hard scatter (no recon- 
structed interaction vertex or fewer than three recon- 
structed tracks), or from events with identified cosmic 
muons. The misidentified jets template is extracted from 
the fake photon sample, which fulfills exactly the same 
requirements as the photon sample except that the pho- 
ton track isolation requirement is inverted. This sample 
is dominated by misidentified jets. Finally, the e/7 tem- 
plate is obtained from a data sample of isolated electrons. 

The total number of background events from misiden- 
tified jets (-/V m isid) can be predicted from the fake photon 
sample based on the rates at which jets, passing all other 
photon identification criteria, fail or pass the track iso- 
lation requirement. To measure those rates we use an 
EM plus jet sample, where the EM object passes all pho- 
ton identification requirements except the track isolation, 
and where the jet approximately balances the EM object 
in the transverse plane. We first determine the number 
of events (N\) in the sample that fail the track isolation 
requirement. We then fit the DCA distribution of the 
events that pass the track isolation to a linear sum of the 
e /7 and misidentified jets templates in order to extract 
the number of misidentified jets (N2) passing the track 
isolation. iV m i s id is then equal to the number of events 
in the fake photon sample multiplied by N2/N1. We fit 
the DCA distribution in the photon sample to a linear 
sum of the three templates, fixing the contribution of 
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FIG. 1: DCA distribution for the selected events in data 
(points with statistical uncertainties). The different his- 
tograms represent the estimated background composition 
from the template fit to this distribution. The inset figure 
compares the individual template shapes. 



misidentified jets as described above, and determine the 
e/7 and non-collision contributions. The result of the fit 
is illustrated in Fig. Q] Most of the signal photons have 
DCA less than 4 cm, therefore we limit our analysis to 
this particular window, which contains 29 data events. 

The only physics background to the j+$ T final state is 
the process Z+7 — > vT> + ~f. This irreducible contribution 
is estimated from a sample of Monte Carlo (MC) events 
generated with pythia u sing CTEQ6L1 parton distri- 
bution functions (PDFs) [13j. The main instrumental 
background arises from W — > ev decays, where the elec- 
tron, due to tracking inefficiency or hard bremsstrahlung, 
is misidentified as a photon. This contribution is esti- 
mated from data using a sample of isolated electrons. 
The same requirements as for the photon sample are im- 
posed, and the remaining number of events is scaled by 
(1 ^ etrk)Atrk) where etrk is the track reconstruction ef- 



ficiency of (98.6 ± 0.1)% [14j. A smaller instrumental 
contribution to the background is expected from W + 7 
production where the charged lepton in a leptonic W 
boson decay is not detected. The kinematics of this con- 
tribution is obtained from VF(+jets) — > lepton + 2/(+jets) 
MC samples generated with pythia, while the cross sec- 
tion is taken from the MC generator based on [l5| , which 
predicts all contributions (initial state radiation, trilin- 
ear gauge boson vertex, and final state radiation) to 
the full process. We generate signal events [16[ with 
M D = 1.5 TeV for n = 2, 3, 4, 5, 6, 7 and 8. For different 
values of Mo, the cross section scales as 1/M^ + , leav- 
ing the kinematic spectra unaffected for a fixed number 
of extra dimensions. 

All MC events are passed through a detector simu- 
lation based on the geant 17] package, and processed 
using the same reconstruction software as for the data. 
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TABLE I: Data and estimated backgrounds. 



Background 


Number of expected events 


7i A- 'v - — > ijIj -4- 'Y 


12.1 ± 1.3 


W -> ev 


3.8 ±0.3 


Non-collision 


2.8 ± 1.4 


Misidentified jets 


2.2 ± 1.5 


W + -J 


1.5 ±0.2 


Total Background 


22.4 ±2.5 


Data 


29 



Additionally, wc apply scale factors, with values ranging 
from 94% to 98%, to account for the differences between 
the efficiency determinations from data and simulation. 

The main sources of systematic uncertainty are the 
uncertainty in the photon identification efficiency (5%), 
the uncertainty in the total integrated luminosity (6.1%), 
and the uncertainty in the signal acceptance from the 
PDFs (4%). 

For the SM backgrounds estimated from MC, the 
quoted uncertainties include the uncertainty in the the- 
oretical cross section, which is dominated by the uncer- 
tainty in the next-to-leading-order K factors (7%). For 
the range of pt in question and for the selection require- 
ments used in this analysis, the K factors vary around 
unity within this uncertainty margin [HEH. The uncer- 
tainty in the width of the e/7 sample DC A template re- 
sults in an additional systematic uncertainty of 0.4 events 
in the non-collision background estimate. 
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FIG. 2: Photon pr distribution for the final candidate events 
(data points show statistical uncertainties), after all the se- 
lection requirements. The LED signal is stacked on top of SM 
backgrounds. 

The final numbers of events for data and backgrounds 
are given in Table |TJ Fig. [2] shows the photon px distri- 
bution, with the SM backgrounds stacked on top of each 
other. Data and the SM expectation agree, so we pro- 
ceed to set lower limits for the fundamental Planck scale 
Mb- We employ the modified frequentist approach [l9[ 



to set limits on the production cross section for the sig- 
nal. This method is based on a log-likelihood ratio test 
statistic and uses the binned photon px distribution. As- 
suming the leading-order theoretical cross section for the 
signal, we derive the following lower limits on Md at 
the 95% C.L.: M D > 884, 864, 836, 820, 797, 797 and 
778 GeV for rt = 2, 3, 4, 5, 6, 7 and 8 extra dimensions, 
respectively. Table [XT] and Fig. [3] summarize the results 
for the limit calculations. 

To conclude, we have conducted a search for LED in 
the 7±$ T channel, finding no evidence for their presence. 
We have set limits on the fundamental Planck scale, sig- 
nificantly improving results of previous searches. 



TABLE II: Summary of limit calculations. 


n 


Signal 


Observed (expected) Observed (expected) 




efficiency 


cross section 


Md lower 








limit (fb) 


limit (GeV) 


2 


0.49 ± 0.04 




884 (921) 


3 


0.4 


3 ±0.04 


24.5 (22.7) 


864 (877) 


4 


0.47 ± 0.04 


25.0 (22.8) 


836 (848) 


5 


0.43 ±0.04 


25.0 (24.8) 


820 (821) 


6 


0.50 ±0.05 


25.4 (22.3) 


797 (810) 


7 


0.49 ±0.04 


24.0 (23.1) 


797 (801) 


8 


0.52 ±0.05 


24.2 (21.9) 


778 (786) 
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FIG. 3: Expected and observed lower limits on Md for LED 
in the 7±$ T final state. CDF limits with 87 pb _1 of data 
and the LEP combined limits are also shown. 

We thank Stephen Mrenna for his help with generat- 
ing MC signal events, the staffs at Fermilab and collab- 
orating institutions, and acknowledge support from the 
DOE and NSF (USA); CEA and CNRS/IN2P3 (France); 
FASI, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST 
(India); Colciencias (Colombia); CONACyT (Mexico); 
KRF and KOSEF (Korea); CONICET and UBACyT 
(Argentina); FOM (The Netherlands); STFC (United 



7 



Kingdom); MSMT and GACR (Czech Republic); CRC 
Program, CFI, NSERC and WestGrid Project (Canada); 
BMBF and DFG (Germany); SFI (Ireland); The Swedish 
Research Council (Sweden); CAS and CNSF (China); 
and the Alexander von Humboldt Foundation. 



[aj 
W 
M 
[d] 



[f] 
[t] 

[1] 

[2] 
[3] 
[4] 

[5] 



Visitor from Augustana College, Sioux Falls, SD, USA. 
Visitor from The University of Liverpool, Liverpool, UK. 
Visitor from ICN-UNAM, Mexico City, Mexico. 
Visitor from II. Physikalisches Institut, Georg- August- 
University, Gottingen, Germany. 

Visitor from Helsinki Institute of Physics, Helsinki, Fin- 
land. 

Visitor from Universitat Zurich, Zurich, Switzerland. 
Deceased. 

N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Lett. 
B 429, 263 (1998). 

G. Guidice, R. Rattazzi, and J. Wells, Nucl. Phys. B544, 
3 (1999). 

D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 
89, 281801 (2002). 

B. Abbott et al. (DO Collaboration), Phys. Rev. Lett. 
86, 1156 (2001); ibid 90, 251802 (2003); ibid 95, 161602 
(2005). 

D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 



[6 

[7 
[8 

[9 

[10 

[11 
[12 

[13 

[14 

[15 
[16 

[IT 

[is; 

[19 



92, 121802 (2004). 
LEP Exotica 



Working Group, 



URL: 



http : //lepexot ica . web . cern . ch/LEPEX0TICA/notes/2004-03/ed_n 

and references therein. 

V.M. Abazov et al. (DO Collaboration), Nucl. Instrum. 
Methods Phys. Res. A 565, 463 (2006). 
Pseudorapidity is defined as — ln[tan(0/2)], where 8 is 
the angle between the particle and the proton beam di- 
rection. 

S. Abachi et al. (DO Collaboration), Nucl. Instrum. Meth- 
ods Phys. Res. A 338, 185 (1994). 
T. Andeen et al, FERMILAB-TM-2365 (2007). 
G. C. Blazey et al, |arXiv:hep-ex/00050i2| (2000). 
T. Sjostrand et ah, Comput. Phys. Commun. 135, 238 
(2001). 

J. Pumplin et al., JHEP 0207, 012 (2002); D. Stump et 
al, JHEP 0310, 046 (2003). 

V.M. Abazov et al (DO Collaboration), Phys. Lett. B 
659, 856 (2008). 

U. Baur and E.L. Berger, Phys. Rev. D 41, 1476 (1990). 
Stephen Mrenna, private communication. 
R. Brun and F. Carminati, CERN Program Library Long 
Writeup W5013, 1994. 

U. Baur, T. Han, and J. Ohnemus, Phys. Rev D 57, 2823 
(1998). 

W. Fisher, FERMILAB-TM-2386-E (2007); T. Junk, 
Nucl. Instrum. Methods Phys. Res. A 434, 435 (1999); 
A. Read, CERN 2000-005 (2000). 



